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ABSTRACT 


The ability to predict torsional instability in the 
early stages of design can Rave important consequences 
on the design of both coriventional and high performance 
ships. This thesis develops fast approximate methods of 
tersional buckling analysis for particular application 
im the concept, feasibility, and preliminary stages of 
ship design. 


Two simplified mocels of stiffeners commonly used in 
ship construction were presented. The first was an 
ideal I section stifferer. The second, a flat bar 
stiffener, which was included in this analysis because 
of its favorable producticnal properties and wide 
useage. These models, with initial imperfections, were 
subjected to axial compressive end loading and the 
resultant behavior analyzed. The approach of the 
analysis of this thesis was an application of the 
energy method to determine the critical buckling stress 
instead of the more commonly used equilibrium aporoach. 
Both beam theory and thin plate theory were used in 
energy-work relationships in these derivations. Then 
the first yield lead was determined utilizing the 
Perry—-Robertsoan approach. Integrated into all phases of 
this exploration was the concept of the geometric 
imperfections of the stiffeners. 
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NOMENCLATURE 


ra Nondimensional parameter 
a Length of stiffener between transverse supports 
b Uniform stiffener spacing 


b Plating effective width 


B Angle of rotation of the stiffener about the toe 


B Initial imperfection of the stifferer in terms of a 
rotation angle about the toe 


B Amplitude of the sinusoidal function for B 


B Anplitude of the sinusoidal function for & 
C Rotational spring constant (moment/length) of the 
supporting plating 


5 Longitudinal warping constant of stiffener 


D Flexural rigidity of plating 


D Flexural rigidity of web plating of stiffener 

d Depth of flat bar stiffener; overall depth of Tee 
stiffener 

d Depth of stiffener to midthickness of flange 


d Depth of web 


S Young’s Modulus of the material 





2 Membrane strain 


f Width of stiffener flange 


G Shear Modulus of the material 


h Height of neutral axis of plate-stiffener 
combination from midplane of plating 


I Effective vertical moment of inertia of stiffener 
and associated effective width of plating 


I Polar moment of inertia of stiffener about toe 
p 

i Vertical moment of inertia of stiffener alone about 
F toe 

I Moment af inertia of stiffener about web plane 


J St. Venant’s torsion constant for stiffener 
= Reduced slenderness ratio 


M(x) Moment developed in the stiffener by the total 
torque 


m Mode number 

p Axial end load 

S Height of stiffener shear center above toe 

§ Elastic axial torsional buckling stress in 


CR 
stiffener (critical stress) 


§ Elastic axial torsional buckling stress limit of 
el. 
stiffener (classical stress) 


§ Axial stress in stiffener 
e 





5 
x 


(S$ 


Axial stress component 


Maximum compressive stress developed by 


x MAXT 


the torque 


Tensile yield stress of material 


Amplitude of the sinsoidal torque function 
Total torque developed in the stiffener structure 
by the loading 


St. Venant’s torque developed in the stiffener 


Warping torque developed in the stiffener 


Plate thickness 


Stiffener flange thickness 


Stiffener web thickness 


Displacement in the longitudinal (x) direction 

Total strain energy of the structure 

Sideways flexure (y) of stiffener 

Initial imperfection of stiffener in horizontal (y) 
direction 

Total work of the applied force 


Vertical flexure (z) of stiffener 





Ww Initial imperfection of the stiffener in the 
vertical (2) direction 


2 Height of stiffener centroid above toe 


VY Poisson’s Ratio 


Note: x,y or z subscripts indicate partial derivatives 


with respect to those coordinates . These are used 
* - oe 


particularly with v, v, w, w, 8B, B. 
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INTRODUCTION 


The subject of this thesis is t3 develop a 
Simplified analytical method of predicting torsional 
buckling of stiffeners used in ship structures. 
Torsional buckling or instability is the phenomena by 
which a column, at a certain axial load, fails suddenly 
in a combined mode af twist and lateral bending of the 
cross section (reference 7). Figure (@-1la) (from 
reference 6) illustrates the failure mode referred to 
here as torsional buckling - in this case for a flat 
bar stiffener. Figure (W-ib) illustrates the failure 
mode of the tripping phenomenon which is not covered in 
this analysis. Figure (8-2) (reference 6) is the 
experimental load deflection curve for the failure of a 
Flat bar stiffener under an axial load. Point A of the 
solid (experimental) curve of fiqure (8-2) is the 
ultimate strength load. As can be seen from figure (a- 
2), the curve proceeds in a downward fashion - the bar 
looses stiffness and fails from this point. 

The ability to predict torsional instability 
in the early stages of design can have important 
consequences on the design of both conventional = and 
high performance ships. It is the intention of this 


thesis to develop fast approximate methods of torsional 


up 





Figure (@-ia) Torsional Buckling 


Figure (@-1b) Tripping 


Figure (@-1) Characterization of Torsional Buckling 
and Tripping in a Flat Bar Stiffener 
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buckling analysis for particular application in the 
concept, feasibility, and preliminary stages of ship 
design. 

To this end, two simplified models of stiffeners 
commonly used in ship construction will be presented. 
The first is an ideal I section stiffener. The second, 
a flat bar stiffener, which is ineluded in this 
analysis because of its favorable product ional 
properties and wide useage. These models, with initial 
imperfections, will be subjected to axial compressive 
end loading and the resultant behavior analyzed. The 
approach of the analysis of this thesis will be an 
application of energy methods to determine the critical 
buckling stress instead of the more commonly used 
equilibrium approach. Both beam theory and thin plate 
theory are used in energy-work relationships in these 
derivatians. Then the first yield load is determined 
utilizing the Perry-Robertson approach as described in 
references (9,14). Integrated into all phases of this 
exploration will be the concept of the geometric 
imperfections of the stiffeners. 

It is the intention of this thesis, using the above 
mentioned approach, to lay groundwork which would have 
the potential for future expansion to include mode 


interaction - to be able to easily analyze combined 
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loadings etc.. 

Particular emphasis is placed here upon § inclusion 
of the geometrical initial deformities of the 
stiffeners which occur during the manufacture and 
through day to day usage of the ship. It 15 well known 
that these geometric imperfections can drastically 
reduce the strength of the structural member. But 
current design standards (references 14,11 etc) have 
not tried to quantify these values. Some design codes 
(example: reference 18) have specified tolerernce limits 
of various geometric parameters which must be met 
before the design formulae are considered valid. 
Granted, in the concept, feasibility ard preliminary 
Benes of ship design where the ship only exists. on 
paper, the designer would have no knowledge of the end 
resulting imperfections of the actual ship. Hopefully, 
in the future, with the advent of better production 
control in shipyards and the incursion of the computer 
into ship production facilities, this information will 
be collected and statistically analyzed. 

In addition, once the ship is manufactured, these 
formulae will provide a quick and easy evaluation of 
the actual critical buckling stress in terms of 
physical parameters. Thus, this analysis has made a 


special effort tec define the imperfections in easily 


1S 





measured terms. This analysis measures the horizontal 
or vertical movement of the shear center of the 
stiffener at its point of greatest deflection and 
assumes a simple sinusoidal distribution. Knowledce of 
the web depth easily converts the imperfection into a 
rotational angle. Here, the stiffener is considered to 
rotate about its toe or base. In other words, the base 
plating is considered much larger and stiffer. Figure 
(@-3) illustrates this simplified concept for a Tee 


stiffener. 
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Finure (9-3) Tee Stiffener 


Currently the method of handling the stiffener 
failure problem is incorporated into the use of factors 
of safety. But, with recent attempts to further 


optimize ship structures - in particular to minimize 
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the amount of weight of a ship that must be devoted tos 
its structure - better, yet reliable, methads must be 
devised. Remember, with a ship, less weight in the 
structure equates to more weight for the paylaad be it 
weapons for amilitary ship or cargo for a civilian 
ship. 

This thesis concludes with a comparison between the 
formulae derived from the simplified models presented 
here and publisned formulae, published finite element 


analysis, and experimentally derived results. 
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1. DESCRIPTION OF MODELS 


Torsional buckling or instability is characterized 
by a twisting of the stiffener about its line of 
attachment to the plating. This deformation pattern 
involves both sideways and vertical flexure (v,w) = and 
rotation (B) of the stiffener as shown in figure (1) as 


described in reference (3). 





Figure (1-1la) Coordinate System 
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ae eve SD 


Y 





Figure (i-1b) Undeformed Figure (1i-1c) Deformed 


Figure (1-1) Characterization of Torsional Buckling 
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1,1 Model 1. 

The first model considered is a simplification of 
an ideal I-section. It consists of an I-beam with zero 
tersional rigidity. The web of the I-beam transmits 


only shear forces. In effect, it looks like two urequal 


flanges: 


eae P77 TT | 
A = 


I~beam Stress Distribution 


Figure (1-2) Model I 


This model has been considered in reference (4). 
However. reference (4) used an equilibrium approach and 
initiel imperfections were ceansidered as ari 
eccentricity of the axial load. The resulting equations 
of the critical stress were quite complicated. 

Figure (1-3) illustrates the geometric parameters 


of this model. 
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Figure (1-3) - Geometric Torsional Buckling 


Darameters for Macel <= 


I - Moment of inertia about the web plane (top flange 


Zz 
only) 


3 
ioe (tf J/1le 
rd fF OW 
s - Height of shear center above toe (origin) (top 
flange only) 


Ss =d 
Cc 


C - Longitudinal warping constant about shear center 


~ 
of stiffener alone (top flange only) 


C = @ 
w 
I - Vertical moment of inertia about toe (top flange 
< 
only) 
= e 
ro~-w ft (d + (i/ledt >) 
t wfc f 
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I - Polar moment of inertia about toe 


J - St. Venant’s torsion constant (tap flange only) 
3 
J = (Ff € d/3 
wf 


a4 









1of Siege ai 








~— 


“skew? ots WV ‘ 
¢ . 
Bt # 9b 
> 


1.2 Model II. 


The second model considered is the simple flat bar 


stiffener illustrated in figure (1-4). 





Figure (1-4) Model II - Flat bar stiffener. 


Figure (1-5) illustrates the geometric parameters 


of this model. 
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Figure (1-5) - Geometric Torsional Buckling 
Parameters for Model Il 


I - Moment of inertia about the wed plane 
Z 
(web only) 
3 
I = (d t )/i1e 
Z Ww W 


s ~- Height of shear center above toe (origin) 
(web only) 


s =d /e2 
Ww 


C - Longitudinal warping constant about shear 


center of stiffener alone 


I - Vertical moment of inertia about toe (web only) 
t 3 
Y = (t d d/3 
t Ww Ww 


I - Polar moment of inertia about toe 


2s 





z —- Height of centroid above toe 


J - St. Venant’s torsion constant 
3 
eee Guste) / 3 
WoW 


D - Flexural rigidity of the web plate 


= 
y 


) 


12(1- y 
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2. MODEL I —- TORSIONAL BUCKLING UNDER AXIAL LOADS 


Since Model [ basically consists of twa unequal 
flanges with the top flange rotating about the bottom 
flange, the torsional rigidity 1s considered only for 
the top flange. Thus the strain energy equation 
consists of three parts (reference 1). The first term 
represents sideways bending. The second term represents 
longitudinal warping. The third term represents the 
torsional rigidity of the top flange alone. The strain 


energy for a length a is: 


a 2 2 c 
V = {/e f (EI v + EC B + GJB )dx 
Q Z XxX WwW XX x 


Also from geometrical properties we have (see figure 


1-3) 3 


Thus: 


a 22 2 2 
V= 1/2 f {EI s B + EC B + GJB ) dx 
YW 2 xx WwW XxX x 


< 
il 


a 2 2 = 
1/2 f EC(I gs + C€C )B + GJB J dx 
B Z Ww xX ™ 


2s 





The following choice of B ig made, in order to 


meet certain Boundary conditions (ie B = 8 at x=, 
xX 
x=a ) and to keep B simple. 


B= 8 sin (m/x/a) 
Oo 


Thus for Model I: 
2 e 2 2 


V= 1/4 a B (mit/a) CECI s + C )(m7%/a) + GI] 
o = W 


fa 
Mm 


ae > ae = 
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Before development of the virtual work equation, a 
note should be said about the sign convention used 
here. Compression is held ts be a positive stress”) and 
tension as negative. The following figures illustrate 


the sign convention used here. 
: 


, —> 
q 


Figure (2-1) Sign convention for inplane axial and 
shear stress components. 
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Figure (2-2) Application of sign convention for 
case of an axial load. 
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Model I is not assumed to be perfect but rather to 


have initial imperfections. These initial imperfections 
% +* t 
v ,w and & and the additional deformations v,w and B 


are illustrated in figure (2-3). 


va 
epee 


amet) 
4 —T 


@ - —lL 
-. a 
s\n P 
SI B 
\I 
Y IN 
megure (2-3) Illustration of initial imperfections 


V¥, wt, B¥ and additional deformations v,w,B. 


The virtual work for the case of an inplane axial 


load becomes (reference 1): 
W = Sf S uly,2) dy dz 
A x 


Thus using figure (1-1), it can be shown that the 
strain of the centroidal axis in the bent 


configuration, &@ (membrane strain), is (reference 1): 
x 


2 2 
@ wm-u + i/e(v ) + 1/2 (w ) 
x x x x 


ey 





for a perfect I section. With imperfections, we 


substitute: 


* 2 * 2 
& = -yu + L/e@tv tv ) + L/2lw tw) 
x x x »™ x x 


using the inextensibility assumptions, set: 


& = Q 
x 
we haves 
*¥ 2 * 2 
u = l/2t(v tv ) + 1/2tw tw ) 
x x x 


a * 2 * 2 
u = i/e i C(v tv ) + (Cw tw) J dx 
Q x xX x x 


where u = uly,2). 


From geometrical relationships (see figure (2-3)) 


v = zB w = -yB 

# Ps Es “ 
vy = 2B w = -yB 

(oa 

v = zB w 2 -yB 

x x x x 

+ * it *% 
v = zB w = -yB 

x x x x 


Substituting and integrating we get: 


30 





B x x xX x 
oe = #2 cote 
+y B + ey BSB + y B ) dx 
x x x x 
a e 2 *e2 
u(y,z) = 1/2 f (z + y )(B + Bd) dx 
Q Xx x 


Again using the simple choice for B and B 


B= B sin(m? x/a) 


oO 
x * 
B = 8 sin(m7 x/a) 
oO 
Thus s:- 
a 2 e 
ufy,2) = 1/2 4, (z ty )C(m7/a)B cos(m 7 x/a) 
Q o 
% = 
+ ((m7/a)B cos(m]¥x/a)] dx 
oO 
a ee * 2 
u(y,z) = 1/2 J (2 +y )(B +B) 
Q o 860 
2 
* Eim7/a)cos(m7x/a)] dx 
2 2 * 2 2 
uty,z) = 1/2(z +y )(B +B) (mi? /a) 
o Oo 
5 
* Cx/2 + (sintl@m7 x/a)/ (4m /a)] 
@ 
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2 2 coe 2 
u(y,z) = 1/4(z +y )(B +B ) (m@/a) a 
o oO 


The work equation now becomes: 


are * 2 = 
W= Sf ® €1/4(z +y )(B +B ) (m77*/a) al dy dz 
A ©é 


o 8 6C 
For the current situation of erd loading, $ = 8 
x e 
coc 2 aoe 
W= 8 C1/4(B +B ) (m/7/a) al ff tz +y ) dy dz 
e o 9 mi 


But note that: 


= = 
ak (zZeetey ) dy dz = I 
R 


Thuss 


* 2 2 
W= 1/4 &® (B +B ) (miF/a) a I 
@ o o p 
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Applying the calculus of variations with respect 


to B to the work and strain energy equations yields 
oO 
(reference 35): 


* 2 2 
m= 172 (8 4B ) Gm /a) at SE 
e o 0 p oO 

2 2 2 

SV =1/2 aB (miF/a) CE(\1 s + C )tmit/a) 

o Zz A 

+ GIISB 

OQ 


Applying the principle of minimum potential energy: 
Wee bv ~Sw=2 4 
we get: 


Sw=dv 


Solving for & : 


2 
2 _ 2 
CECI gs + C imu /ad + GTI B 
2 Ww © 
& ee 
e * 
I B + 8B 
p o Oo 
(2-2) 


It can be seen by inspection that the lowest 
buckling stress occurs for one wave, m= 1, since m 


must be an integer. Thus: 
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§ =6§ 


g CR 
As it can easily be seen, as B becomes large $ 
| o CR 
approaches the limit of: 
2 2 
CE(I s +C )(m/¥/a) + GJ] 
z W 
3 I ee eae ee eee C2—3) 
CL 
I 
p 
Figure (2-4) illustrates the behavior of §& with 
t CR 
increasing B +8. 
oO ° 


Further simplification and rearranging of equations 


(2-2) and (2-3) provides: 


re eS (2-4) 
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Figure (2-4) 8 vs. (B + 8B) 
CR O oO 


33 





e.9 Determination of the First Yield Load. 

To determine the first yield load, the total torque 
developed in the torsional buckling process must be 
first evaluated. This torque will consist of two parts. 
The first term consists of the well known St. Venant’s 
Torque (reference i). The second term represents the 


warping torque (reference @2). 


From reference (1): 


where: 


T = GJB 
SV x 


and from reference (2): 


T = - ET dew 
W 2) = Xxx 


and from geometry; 


resulting in: 
e 
feo= — EL gd B 
* ZC Xxx 


Thus the total torque is: 
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2 
T= GJB - EldB (2-5) 
x <= C XXX 


Figure (2-3) illustrates the above equations. 


Again a simple choice of B is made, and with m= 1: 


B= B sin(? x/a) (2-6) 
oO 


Equation (2-3) now becomes: 


T = GI(7/a)B cos(7?x/a) 
O 


2 3 
+ EI d (7*/a) B cos(7?x/a) 
Zc oO 


2 2 


T = (7 /a)B cos(*x/a) (GJ + EI d (7 /a) ] 
O ze 


Choosing a simple sinusoidal shape for T aliows: 


T= T cos(7 x/a) 
oO 
Thuss 
2 2 
T = (¥/a)B (GI + EI d (sa) 3 (2-7) 
oO oO 2c 


To eliminate B , equation (2-4) is applied to equation 


o 
(2-7). 
* 
B 
Oo 
2 2 
T = (63 + EI d (#7 /a) 1(%/a) ------------- 
o - ¢« 
(§ /8 —- 1) 
El eR 


(2-8) 
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From sectional considerations, an evaluation of the 


Maximum stress can be made. 


Figure (2-5) Sectional illustration. 


The moment developed by the torque in terms of 
deflections is (reference 2): 


Mex) = - EI w. 
2 xx 


So the maximum compressive stress developed in the 


Flange due to the torque is: 


- EI dg B f 
ZC XxX Ww 
(& ) Se a re ee co ee co ee ree ee ae ee me ne 
x MAXT 
el 
z 
(S ) =- Ed B f /2 
x MAXT Cc XxX w 


Application of the above definition for B (equation 


2-4) leaves: 
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Ea f 


cow e 
(S ) = ---------- (7 7a) Besime xy a) 
x MAXT Oo 
2 
A choice of x iS made such that (8S ) 1S a maximums: 
x MAXT 
Ed f 
cow 2 
(2) ) = ---------- (W7/a) B (2-9) 
x MAXT Oo 
2 


Equations (2-7) and (2-9) are combined to eliminate B. 


oO 
Ed f T (7/a) 
cwo 
(S ) = see (2-19) 
x MAXT 2 2 
2(GJ +EI d (7 /a) 3] 
Ze 
Using equation (2-8) leads to: sa 
2 re 
Ed f (i /a) B 
cow oO 
(§ ) Senn nnn nn oe (2-11) 
x MAXT 
2 cS /§ - iJ 
CL CR 
% 
Note that when B = 8 that this equation leads to an 
o 
indeterminate form of 8 since & = 8 , 
4 CL CR 


The total maximum stress can now be described by: 


(& ) = (§ ) + D/A 
x TOTAL x MAXT 


where 
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P/A = § 
CR 


Thus the yield stress has been shown to be: 


S§ = (§ ) + & Cea) 
y x MAXT CR 


Using equation (2-11) in equation (2-12): 


2 << 
Ed f (7 /a) B 
Cow 
Bm mmr wr eae “re = 
y CR 
= & /% ae 
CL CR 
C213) 
2 
Solving for & and dividing by & to place into 
CR y 
nondimensional form yields: 
= c * 
es & Ed f (isa) B & =f 
CR CL. cow o CR oe 
-—-—— - Ll tom FO enn nr nr nr rrr -—— ---|/= @ 
§ § 2s & & 
y y 7 y y 
(2-14) 


Choosing the smallest root and simplifying the equation 


by introducing a slenderness ratio L: 


ve=12) 
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2a 
7 


‘et cf euocig ose €2689 





and letting: 





= * 
(7 /a) Ed f B 
C W O 
A= e3rree------- (2-16) 
2s - 
y 
Then ¢: 
= « 
§ (1 + L/L +€A) —- (1 + {/L + A) - 4/L 
CR 
& 2 
y 


teo=17) 

A plot of § /8§ vs L for various values af A is 
shown in egret Figure (2-6) was compared with 
various design criteria from reference (18) and it was 
found that equation (2-17) was similar to design 
equations of reference (19). However, the reduced 
Slenderness ratio, L, of reference (18) is based upon 
simple Euler buckling and the geometric parameter, A, 
is of course different. The colunm curves of reference 


(12) were quite similar in appearance and range of 


values to that found in figure (2-6). 
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Figure (2-6) 
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3. MODEL II - TORSIONAL BUCKLING UNDER AXIAL LOADS 


Since Model II is the very simple case of a flat 
bar stiffener, thin plate theary can be easily utilized 
in the development of the strain energy equation. Much 
of this section has been patterned after a similar 
analysis presented in reference (3). 

For a web plate of uniform thickness and a length a 


the strain energy is repensented by (reference 5): 


a (x) c 
V= i/7e D 4 Civ + vy ) 
o xx Z2 


W oo 
2 
~2(1-V)Ctv v-v J] dx dz 
XX 22 Ke 
a 2 
+ 1/2 [ CB dx 
oO 


The second term represents the rotation of the 
supporting plate structure modeled here as an elastic 
spring. D is the flexural rigidity of the web plate. 
From Pe ical considerations: 

Vv = 2B 
Making a Simple choice for B; 

B= 8 sin (m7? x/a) 

o 


substituting into the expression for the strain energy 


and integrating over a constant depth d results in: 
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2 2 2 Zz 
3S) tel Cymer eee) Sl Ys ale 
Ww © 


=~ 


e 
+ CB a/4 (3-1) 
Oo 
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As for Model 1, Model II is not assumed to be 

perfect but rather to have initial imperfections. These 
~ * Ps 

initial imperfections v, w, 5 and the additional 


deformations v, w, andB are illustrated in figure 


moni). 





Figure (3-1) Illustration of the initial imperfections 
V¥, wt, B® and additional deformations v,w,B. 


For the case of an inplane axial load, the work 


equation is : 


W = Jf & uly,z) dy dz 
A xX 


From figure (3-1), it can be shown that :; 
a 2 
&@ =z2-—-u + 1/2(v ) + 1/2w ) 
x x x x 
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for a perfect web section. With imperfections, and 
using the inextensibility assumption with regard to the 


length of the stiffener setting ® = 9, the result is : 


x 
* 2 * 2 
uo = 1/2 (v +v) + 1/e@lw +w ) 
x x x x x 
a ae * 2 
u = i/e f Civ +v) + 1/2e@€w + w) J dx 
o x x x x 
(3-2) 
where u = uly, 2Z) 
As shown in figure (S=1)5 the geometrical 


relationships (equations 2-1) are used to simplify 


equation (3-2): 


a e 2 * 
ufy,z) = 1/2 { (z + y ){B +B) dx 
G 


“t 
A simple choice is made for B and B. 


B = 8 sin (m/x/a) 


oO 
% + 
B = 8B sin (mifx/a) 
oO 
Thus: 
2 2 + 2 
u(y,2) = 1/4(z + y )(B + B )d(mi/a) a 
oO O 


The work equation now becomes: 
= 2 * 2 


W = Sf S&S C1/4(z + y )(B + Bd) #/a) al dy dz 
A x o G 
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* 2 2 
W= 1/4 8 (B +B) (m@/a) ail 
x oO S p 


Applying the calculus of variations with respect to 


B to the work and strain energy equations yields: 
o 


* 2 2 
ois 1/28 (Bo +B) (maa) at SB 
x oO o Dp oO 


2 2 2 
{D ad im”“/a) B Cd (m/a) 
W o 


Sv 


2 
+ 6(1 -Y) 1/6 + CB a/2 $B 
oO oO 


Applying the principle of minimum potential energy: 
Tl=Sv- Sw=a@ 
SW = dV 


Solving for § =8& 3: 


e CR 
an 2 
Ddim™@/a) Cd (m*/a) + 6(1-V)I/3 + C B 
we o 
4 ms eee ee a i i i ee ee ee 
CR = 2 a 
(m4 /a) I Bots 
=) o o 
(3-3) 
For the case of zero restraint against reotatican; 
C= @: 
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ae 2 


D d(m/?/a) Cd (m*/a) + 6(1-Y)]} EB 
Ww o 
ee ee a a ca ee a a a a ee ce cs ee ee ee ee ee ee en ee ee ee 
CR + 
a B +8 
p oO Q 
(3-4) 
It is easily seen that the lowest §& for this case 
CR 


results from the simplest of mode shapes, m=1. 

Wher the rotational restraint 1s not equal to zero; 
the lowest value of § depends upon the degree of 
rotational restraint: - 

3 1/4 
m= (a/iF)(3C/D d ) (3-5) 
W 

With the use of equation (3-5), m may be evaluated 
ayjd rounded up or down to the nearest integer for the 
correct wave shape. 

It is easily seen that as B becomes large ae 

o 


approaches a limit. Thus with (8 ) =S5 ¢ 
CR LIMIT Cit 


a = 
Dd(mi/a) Cd (mi*/a) + 6(1-V)1/3 + C 
Ww 
& Ss a nn en a a = 
et = 
(m 7a) I 
p 
(3-6) 
Figure (3-2) illustrates the behavior of §& with 
+t CR 
increasing B +8. 
oO o 
Further Simplification and rearranging of 
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equations (3-3) and (3-6) provides: 


RM (3-7) 


S 





= 
7 
Bo Bo + Bo 
t 
Figure (3-2) Behavior of §& with B +B 
CR ° oO 


The first yield load will not be determined for 
this case under this method of analysis due to the 
relative complexity of equation (3-3) and due to the 
fact that the Perry-Robertson method will not lend 
itself to further simplification along this train of 
analysis. Since a simplified formula for design 
purposes is the goal of this thesis, the next chapter 


will use a different tactic to solve the probienm. 
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4% MODEL If —- GENERALIZED ANALYSIS 


aS => ae eee oS Se a SS —_ee ae ee oS =e ee Se Se ae Ss ee 


The analysis thus far for the case of a flat bar 
stiffener under end loading has proceeded alorg the 
lines of thin plate thecry up to the development of an 
expression for the critical buckling stress. 

This chapter uses beam theory to develop a simpler 
formulation ina manner analagous to that performed on 
Model I in chapter 2. In fact, this analysis will 
merely be a generalization of the results for Medel iI 
applied to Model II. It must be remembered here that 
Model i consisted of two unequal flanges with the top 
flange deforming or rotating about the bottom flange. 
At no time in the derivation of the formula for the 
critical stress of the tee stiffener is there an 
assumption made about the position of the top flange 
relative to the bottom flange. Thus the same analysis 
with minor changes will be performed for the flat bar 
stiffener. However, in the case of the flat bar 
stiffener of Model II, there is a restraining effect 
provided by the jincture of the flat bar with the 
Plate. This in effect can be modeled as a Spring with a 
rotational constant C. Thus the strain energy equatisn 


would have four terms (instead of the three terms for 





Model ID. These terms represent (in order of 
appearance) sideways bending; longitudinal warping; 
torsion; and rotation of the supporting plate structure 
modeled as an elastic spring. Thus the strain energy 


equation is: 


a = 2 o 2 
Ve= i/2 f CEI v + EC B + GJB + CB J dx 
oO 2 XX W XX x 


Again choose: 


B = B sin(m itx/a) 
© 


Substituting and integrating results in: 


0 


2 2 4 2 
V=1/4a B CECI s +C )(mir/a) + GIt(mim/ad) + C ] 
o a W 


Applying the calculus of variations with respect to 


B yields: 
© 
@ 4 
RUS 1/2 4 B CECI s +C )tmW/a) 
OQ z W 
2 
+ G6Jim WY/a) *cide (4-1) 
© 


The assumptions used in deriving the work equation 
for Model I are the same as for Model II. Thus’ from 
chapter 2 we have: 

* 2 rad 


W= 1/48 (B +8) (mit/a) a I 
2 o o p 


v1 





and 
KW =i/es ( + B)imP/abe al OB 
e O O p oO 


(4-2) 


Using equations (4-1) and (4-2) and applying the 
principle of minimum potential energy and solving for 


the stress yields (with $§ =8 )s 


e CR 
2 4 2 
E(Is +C )(m¥/a) + GJ(miT/a) +C B 
z Ww Oo 
- Ss meee ee 
CR 2 * 
(m7-/a) I B +8 
p o oO 
(4-3) 
Let: 
2 4 2 
E(I s +C )(mV/a) + GItmP/a) +C 
Pa Ww 
& Sn nn nn rn nr (4-4) 
cl 2 
(m@say I 
p 
Using (4-4) in equation (4-3) and solving for 8 
Oo 
yields: 
+ 
B 
° 
B = ----------- (4-5) 
o 
gs /8% - 1 
Se CR 


Comparison of equation (4-3) with equation (e-1) 


v2 





reveals the only difference is that equation (4-3) has 
an additional term in the numerator (specifically the 
spring constant C). 

Since 1t 1s very difficult to quantify the spring 
constant C and (hopefully) since its value is small in 
comparison to other terms in expression (4-3), C will 
be set to zero value and ignored for the remainder of 
this analysis. 

Thus for the case of C = @, equation (4-3) reduces 


to equation (2-2). 


4.4 Determination of the First Yield Load. 

The determination of first yield lead far the 
generalized case proceeds in a similar fashion as for 
the case of Model Il. 

Again the total torque consists of two parts 
(reference 1), St. Venant’s torque and a warping 
torque. However, in this case the warping torque 15 
(reference 2): 

T += -EI dw 

ws 2 xX 


and 


w = dB 
x x 


resulting in: 
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2 
foeoceeel cd B 
W 2 xX 


Thus the total torque is: 
= 
T = GJB - EI dsB 
x z x 


Choosing B and T: 


B= B sin(mT x/a) 
Oo 


T = T cos(m/ x/a) 
Oo 
Leaves: 
2 c 
T = (#/a) CGI + EI d (iF/a) 1 B (4-6) 
oO = 2 


Using equation (4-3) to eliminate B leaves: 


So 
* 
B 
2 2 oO 
fom 7aytes + Fl d (7'/aj) j ---—-——----———— 
oO = ‘ 
(S /& — t) 
CL CR 


(4-7) 
Again sectional considerations can be used to 
evaluate the maximum stress. The maximum stress 


developed in the stiffener due to the torque is now: 


2 
(S ) = -Ed B 
x MAXT x 


Application of the above definition of B and a 


choice of x such that (8 ) iS a maximum leaves: 
x MAXT 
2 2 
(% ) = Ed (77a) B 
x MAXT o 
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Using equation (4-6) to eliminate B: 
oO 
= 
Ed T (W/a) 


(S ) = -~------------------ (4-8) 
x MAXT 2. 2 
[GJ + EI d (W/a) ] 
< 


Using equation (4-7) in equation (4-8) leads to: 


fe 
B 
2 2 oO 
(& ) = Ed (i /a)  ------------- (4-9) 
x MAXT 
(S /& -— 1) 
CL CR 


Again, from sectional considerations, the total 


maximum stress can be described by: 


S&S = (§ ) + §& (4-19) 
y x MAXT CR ; 
Using equation (4-93) in equation (4-123), solving for 
2 
& and dividing by 8 to place into nondimensional 
CR y 
form yields: 
e e oF 
& & Ed (j*/a) B es 
CR El oO CR ee 
am——/| — jl & e-e- + werrrrrrrere-n- --— + {| ----| = @ 
& & s S & 
y y y y y 


(4=11) 
Choosing the smallest root and simplifying the 
equation by the use of the reduced slenderness ratio 


(equation 2-15) but redefining parameter A as: 


oo 





ace 
(iafay Ed B 


po (4-12) 


Results in (equation 2-17): 


v 


2 2 c 2 
S wee i/o + A) = Ce Ee A = 47 


(4-13) 
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ge COMPARISONS WITH OTHER STUDIES. 


Iza 
ju 
IS 
me) 
jr- 
ify 
|ct 
iD 
[su 
IS 
1a 


In chapter 3, the critical buckling 
flat bar stiffener was derived through the 
plate theory: 


2 2 2 
Ddt(mZ7/a) Cd (m7P/ad) + 641-V)I/3 


O 
AI 
2 


a 


(mt/a) I 


- 


In chapter 4, the critical buckling 
flat bar stiffener was found using one dime 
theory: 


2 4 2 


(mi /a) I 


stress of a 


use of thin 


+C 5 
oO 
* 
Be ts 
oO 2 
(3-3) 


stress of a 


nsional beam 





It should be noted here that whens 


2 2 2 
Ddim7#/a) Cd (m7? /a) + 6(1-V)1/3 = 
WwW 
2 4 2 
E(I s + C )m?/a)d) + GIlm v/a) 
z iw 


the agreement between the two methods would be exact. 


For the elastic range of experimentation and with: 


3 
Et 
w 
D = -~------ G = E/2(1+ V) = E/2.6 
W 
12(1-V) 
3 S 
I = (d t )/12 J = (d t 7/3 
2 w WwW WoW 
d=zd Y= 0.3 
w 
The above relationship reduces to: 
= 2 e 2 
d (mi7i/a) 1 d (mi/a) 1 
Ww Ww 
ee ee ee ee eee + ee a —— ee + — 
10. 92 2.6 le 7.8 
(Thin Plate Theory) (Beam Theory) 


From a cursory examination it would appear that the 
Thin Plate Theory results (equation 3-3) will probably 
Give a larger value for the buckling stress for the 
normal range of stiffener dimensions used in ships. [Ts 


test this hypothesis, the dimensions of a flat bar 


~8 





stiffener used in experiment #21 of reference (6) will 
be used. 
where: 


gd = 87.4 ecm 
Ww 


a 732 cm 


For the thin plate theory: 


For the beam theory: 


e 2 
d (m i/a) 1 


Thus, within the elastic range, the thin plate 
theory results would indeed provide higher buckling 
stress. The beam thoery results may thus be toa 
conservative. However, section 3.3 shows good agreement 
between equation (4-11) (first yield load) and 


experimental yield results. 
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Reference (3) performed a similar analysis to that 
performed in this thesis but without the inclusion of 
initial geometric imperfections. 

For the case of a Tee stiffener subjected to an 
axial load, equatian el of reference (3) allows (in the 
notation of this thesis): 


ce 2 2 


GJ + E€(I s + C )(mi¥/a) + Cla/mir) 


Equation (2-2) of this thesis, for a Tee stiffener 
in an analysis which ignores the contribution of the 
web and includes the contributions of initial defects, 


iss: 


2 2 
CECI s + C )imit/a) + GJI B 
- ” O 
& ma a] = A a ee a a oo oe ee = = ee 
CR * 
I B +8 
p o o 
Thus equation (2-2), for the case of an initially 
perfect Tee stiffener reduces to equation 21 of 


reference (3) when the rotational restraint is held to 
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be zero. However, it must be noted that since the 
physical models of the two derivations were slightly 
different - namely that the model used in this analysis 
consisted of two unequal flanges - some of the 
parameters of the equations will have different values 
ore. [ , J etc.). 

For the case of a flat bar stiffener, equation 22 


of reference (3) has: 


22 2 
DdC€im*/a) d + 6(1-Y)] Cla/m iv) 


From chapter 4, equation (4-3): 
e “ = 
Et(Is + C )(m#/ad) + GI(mir/ad + C B 


< W c 


I (m 7F/a) B +B 

p Oo o 
The following is a comparison for a Tee stiffener 
and a flat bar stiffener which was presented in 
reference (3). Figure (5-1) from reference (3) shows 
the geometric dimensions of the Tee and flat bar 
stiffeners. Table (S-1) shows how equations 21 and ee 
of reference (3) and equations (2-2) and (4-3) of this 


thesis compare with a finite element analysis performed 


el 





using a computer program developed and documented 


the University of California (reference 12). 





4 
Flat Bar (CS /E) x 18 C = @ 


CR 


mye | ee e _ 
pmse] 758 | sk = 









4 
Tee (S /E) x 12 C= @ 


CR 
eeefas [| me | = 







Table (5-1) Comparison with Published Formulae 
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FLAT BAR STIFFENER —a © 100 em 


i= 7.90 em 


14.33 em 





4 
AL 
b, = 30.5 em 


TEE STIFFENER —2a © 160 cm 


Figure (5-1) Stiffener Geometries for 
Comparative Solutions 
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The U.S. Navy has in the past used similar 
formula as those presented here for the determination 
Of critical buckling stresses. For the critical 
buckling stress of a Tee stiffener, in the notation of 
this thesis from reference 13 we have: 


2 
GJ + (m7v/a) EC 


where here: 


32 33 
Ffdt /i@ +det 
Ww Ww fF Ww 


O 
1} 


GJ 


3 2 
Geevee/ Soe 7d feu tft 
p W wow w Ff 


ms 
It 


This 1S very similar to those results determined by 
reference (5), with the exception that the abeve 
formula does not account for the ratational restraint 
provided by the junction of the web and plate. Also I 


is much simplified here and C , the torsion-bernding 


Ww = 
constant, equates toa the quantity (I 5 +) 2G.) atch 
Z W 
reference (3) - but also somewhat simplified. 
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Reference (6) provides insight into the torsional 
buckling process through the experimental results 
published there. This report describes model tests on 
the collapse and post failure strength of Tee struts 
Simulating flat bar stiffeners in a stiffered plate. 
The T struts were subjected to an axial Icad up to and 
beyond failure. The geometrical imperfections of the 
models were also recorded. The researchers of reference 
(6) (of Det norske Veritas) are to be lauded for the 
thoroughness and completeness of their report, 
particularly in the treatment of the experimental data. 

Table (5-2) contains the model data pertaining ta 
this analysis for ten experiments of reference (6). 
Figures (3-2 to 5-11) show the load- deflection curves 
for these ten experiments. The solid line represents 
the experimental data. Point A of each of these graphs 
represents the first yield load point as determined by 
equation (4-11). 

Notice that there is very good correlation between 
the experimental first yield load point and the results 


of equation (4-11). 
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# 
Model a A ae ae ws 


SI ee 
E [ons fee | ares [er Ta 
2 See 
fe [ome fenr [ne [eee [ae | a [oe 
BU fac | a0)| aor-es | ox [ew | = 
[more [ee ees [a Leo | 
[oo fe [ue ees [as [oe Ts 
[ote 


Notations 


nue 





a - length of model (cm) 
bh - height of stiffener (cm) 


t -—- thickness of stiffener (cm) 


Ww 
* 
B - angle of rotation (radians) 
o 
2 
& - yield stress (N/mm ) 
y 
= 
§ - first yield load stress (N/mm ) 
CR 
2 
& - ultimate stress (experimental results) (N/mm ) 


u 


Table (5-2) Data of Model Tests 
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Figure (5-2) 
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Figure (5-3) 
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Figure (5-4) 
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Figure (5-5) 
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Figure (5-6) 
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Figure (5-7) Modei Test #265 of Reference (5) 
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Figure (5-8) Model Test #27 of Reference (6) 
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Figure (5-9) Model Test #28 of Reference (6) 
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Figure (5-19) 
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6 CONCLUSIONS 


This thesis utilized a design approach to tackle 
the problem of torsional buckling of stiffeners with a 
good deal of success. Integrated into all phases of 
this analysis was a consideration of the geometric 
imperfections of the stiffener — both Tee and flat bar. 
These structural members were analyzed as simplified 
sections and the Perry-Robertson approach used in other 
areas of engineering was successfully applied here. 

Chapter 3S showed that the derived results of this 
thesis reduced to currently published formulae for the 
case of the perfect stiffener. The contribution of this 
thesis is that the initial imperfections were taken 
into consideration, using a simple design approach, in 
the derivation of the results presented here. 

Also shown in chapter S was a generally) good 
correlation between results of this thesis and axial 
load tests on flat bar stifferers. Comparison of the 
Tee stiffener formula with experimental data, which 
measures and reports initial imperfections, needs to be 
performed. Also, comparisons should be made with other 
tests on flat bar stiffeners in more dimensional ranges 
to determine the range of validity for the results 


derived here. 


a7. 





Further work needs to be done in this area. 
Specifically, determination of the critical buckling 
stress and first yield load due to lateral Iloadings. 
Also, a Simplified method of handling the meade 
interactions of combined lateral and axial loadings is 
needed. In addition, 1t 1S an asset to understand fram 
the work reported in references (6) and (15) that the 
falling path of the load-deflection curve (figure 6-1) 
1s represented by the tripping mechanism which is neat 


covered in this analysis. 
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